404] 135 


404. 


REPRODUCTION OF EULER'S MEMOIR OF 1758 ON THE 
ROTATION OF A SOLID BODY. 


[From the Quarterly Journal of Pure and Applied Mathematics, vol. IX. (1868), 
pp. 361—373.] 


EULERS Memoir “Du mouvement de rotation des corps solides autour d'un axe 
variable,” Mém. de Berlin, 1758, pp. 154—193 (printed in 1765), seems to have been 
written subsequently to the memoir with a similar title in the Berlin Memoirs for 
1760, and to the “Theoria Motus Corporum Solidorum &c." Rostock, 1765, and there 
are contained in the first-mentioned memoir some very interesting results which appear 
to have escaped the notice of later writers on the subject; viz. Euler succeeds in 
integrating the equations of motion without the assistance furnished by the consideration 
of the invariable plane. In reproducing these results I make the following alterations 
in Eulers notation, viz. instead of zx, y, z I write p, q, r; instead of Mo Mb, Mc 
(where M is the mass) I write A, B, C, these quantities denoting the principal 
moments, and in some equations where the omission or insertion of the factor M is 
really immaterial I write A, B, C in the place of a’, b, c?; moreover instead of Euler’s 


DE ad gue aa h9 | 
A HM. C (which denote respectively et y : a i Aza) I write L, M, N; but in 
other respects Euler’s notation is preserved. The equations of motion are 
Adp + (0 — B) qrdt = 0, 
Bdq +(A — C) rpdt = 0, 


Cdr +(B — A) pqdt =0; 
so that putting for shortness 
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and introducing the auxiliary quantity u such that du= pqrdt, we have 


p=A+2Lu, 

g=8+2Mu, 

r=(f +2Nu, 
where A, B, © are constants of integration, and thence 


du 
arm (B + 2Mu)} (€ + 2N 4)” 
where the integral may without loss of generality be taken from u=0; u, and 
consequently p, q, 7, are thus given functions of t; and it is moreover clear that 
A, B, © are the initial values of p, q? 7. We have also if w be the angular 
velocity round the instantaneous axis 
c? =%+8+C4+2(L4+M+N)u. 

Euler then assumes that the position in space of the principal axes is geometrically 
determined as follows, viz. (treating the axes as points on a sphere) it is assumed 
that the distances from a fixed point P of the sphere are respectively l, m, n, and that 


B 


Di 


A K C 


the inclinations of these distances to a fixed arc PQ are respectively X, pu, v. We have 
then the geometrical relations 


cos? | + eos? m + cos?n = 1; 
? 


du ules cos 7 aiii yju COS M COS n 
p sin m sin n? ju sin m sin n’ 
4 cos m cos n cos } 
sin (y — A) == cos (y —A)=—— zed 
( ) sinnsinl ? ( ) sinnsinl ’ 
1 eos n cos | cos m 
sin (A — 1) = = cos(A— u)2 —————; 
( p) sin lsin m ' ( p) sin l sin m’ 
whence also 
e — cos A cos n — sin X cos } cos m 
sin w= a, 
sin ¿sin m 
sin A cos n — cos X cos l cos m 
COS p= = AA 
sin / sin m 
, cos X cos m + sin X cos [cos n 
SP m0 prin Tra v^ 
sin / sin n 
— gin X cos m — cos X cos! cos n 
CRT d. eee r x, 


sin / sinn 
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The geometrical equations connecting the resolved angular velocities p, q, 7 with 
the differentials of l, m, n, X, p, v are 


di sini =dt(q cosn —rcosm), ddsin?l =—dt(q cosm+r cosn), 
dm sin m=dt(r cosl —pcosn), dusin?m=—dt(r cosn + p cos! ), 


dn sin n =dt(pcosm—qcosl), | dv sin?*^n =—dt(pcosl +q cos m). 


Multiplying the equations of motion respectively by cos/, cosm, cos», and adding, 
we obtain an equation which is reducible to the form 


d (Ap cosl + Bq cos m+ Cr cos n) = 0, 
whence integrating 
Ap cos l+ Bg cos m + Cr cos n. =D, 


D being a constant of integration. One other integral equation is necessary for the 
determination of the angles J, m, m The expressions for dl, dm, dn give at once 


p dl sin 14 q dm sin m +r dn sin n = 0. 
Instead of the ares l, m, n, Euler introduces a new variable v, such that 
v = pcos l + q cos m + r cos n; 
by means of the last preceding equation, we find 
dv = dp cos l + dq cos m + dr cos n, 


and then, substituting for dp, dq, dr, their values, 


dv o= (o + M cos m + N cos *) du, 
P q t 


from which the relation between v and u is to be determined. We have 


eos l+ cos m-F  costn=l, 
Ap cos l+ Bq cos m + Cr cos n= D, 
pcos l+ qcosm 4 rcosn=0, 


which give cosl, cosm, cosn in terms of u, v; the resulting formule contain the 
radical 
INN I M*B*rep* + N?C*p*q?) a > (a? of. y Ep 2?) 
y 29v (Ap + Bg! + 07) - v (A*p*+ Beg + Cr)? 
which for shortness is represented by y{(-)}. We then have 


_ Dp (NCg — MBr?) + BCpv (Mr* — Ne?) + LAqr ví: LA 


cos | LA 224 + M Brp + N*C*p? 7; 
citi aie OE ar DAON cct iR ICE 
PA + M? Bp? + Npe 
ONG Ee o is deti Med der E E EL 
LA gr + MB rp? + Npa 
G VE 18 
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and substituting these values in the differential equation 


dv  Lcosl Mcosm | N cosn 


du p q T 


the equation to be integrated becomes 


x (Ag + MBr + Nepe) = LND (Ap? + Be + Cr?) - LMNv (Ap? + Beg? + Cr?) 


+ x (L?Aq?r? + M*Br?p? + N°Cp*q’). 
Now substituting for p, q, r their values, we have 
BAG + eB rp + Npe = PASE + MBECA + MOAB —2LMNu (AA? + BB + GC), 
BPAgr + MBrp? + NCOpr? — LASBG + WBOA + WCAB —2LMNu (AA + BB + GC), 
P+¢Ptr=A +B +€4+2L4M4+N)u, 
Ap + Be + C? 2 (A. +BB + GC, 
Ap + Bg? + C)? = AA SB + CC: 


and writing for shortness 


A + B+ E — E, 
AA + BB+ GC = F, 
NA" BB+ EC*= E, 
PA BE + MBECA + VCAB =H, 
PASC + MBA + WOAB = K, 
where K = EG — F°, substituting these values and observing that 
L+M+ N =- LMN, 
the radical of the formula becomes 
VIC) 24 (K-2LMNGu+ 2D°LMNu — DE +29 Fo — Go), 


and the differential equation becomes 
dv 1 
25 (K -2LMNGu)- LMNSF —LMNG»v-^ pes (H 2LMN Fu) y {(-)}, 


which can be reduced to the form 


Kdv — LMN F Ddu — 2LMN Gudv + LMN Gudu Hdu — 2LMN Fudu 


V {K— DE + 3LMN ($5 — G)u + 2DFo— G9| ^ J (Lu - 3) GMu 43) 2Nu ©) 


www.rcin.org.pl 


404] ON THE ROTATION OF A SOLID BODY. 139 


Euler remarks that as the right-hand side of the equation contains only the 
variable u, the solution will be effected if we can find a function of u, a multiplier 
of the left-hand side; he had elsewhere explained the method of finding such 
multipliers, and applying it to the equation in hand, the multiplier of the left-hand 
or what is 


: , ; , 1 
side, and therefore of the equation itself, is found to be K-3LWNGw 


the same thing QT 


Multiplying by this quantity, the right-hand side may for shortness be represented 
by dU, so that 
dU (H —2LM NF») y (G) du 


~ (K -2LMNGv) y (2Lu +A) (2Mu + $) (2Nu + €) ' 


and U may be considered as a given function of u, or what is the same thing of t. 


As regards the left-hand side, attending to the equation K = EG — FF, the radical 
multiplied into y (G) may be presented under the form 


V[((G — D’) (K — 2LMN Gu) - (Gv — SPyj]; 
and consequently the left-hand side becomes 


(K —2LMNGu) Gdv + LMNG (Gv — DF) du 
(X—LMNGu) (G — $9) (K - 2LMNGu) — (Gv — SF y] ' 


which putting for the moment K —2LMNGu-gp, Gv-DF=q, G — $?—f?, becomes 


pv (F P-P) fp 


pda . qdp the integral of which is sin 1. hence restoring the values of p, q, f, 
the integral is 


obiit Gv — OF 
y (G — D) y (K -2LMNGu)' 
Hence considering the constant of integration as included in U, or writing 


U-G (H —2LMNFu)  (G) du 
+) arco J la tM Ger GT 


we have for the required integral of the differential equation 


Gv- DF 


sin” (G-S5.KKX-3LMNGw) O 
whence also 
Gv — SF feda. 
J(G—-839JK(KX-2LMNGu) ^ 
and 
VIG — 9») (K — 2LMNGu) - (Gv 9 FY]] _ eos U, 


y (G — D?) y (K -2LMNGu)| 
18—2 
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so that the value of the original radical is 


yy V (G - D) v ((K-—22MN Gu) 
/ i( )- X (G) ———— cos U. 
Substituting in the expressions for the cosines of the arcs l, m, n, these values of v 
and the radical; the formule after some reductions become 


_ DAp , BOP (MG — NB) y (G — D?) LAqr y (G — 3») 
ol 270950. ae Nna a HO LN £5. 


DBy , CAq(NA- LG) /(G— 9»)... MBrp y (G — D) 
008 tt a p UE MN o0 d zar pare SEN OS S U^ 


DOr , ABr(LB— MW)J(G—$*) . y, — NCpqv(G- D) 


TENTE GJ(K-3LMNGs — 9^ U * y (Gy y (K — 314 NGuy 9 Y 


where for shortness p, q,r are retained in place of their values y (2Lu + W), y (2Mu +), 
y (2Nu + 6). 


The values of J, m, n being known, that of A could be determined by the 
differential equation 


dt (q cos m + 2 cos n) 
sin? 


NES 


> 


and then the values of m, v would be determined without any further integration; 
but it is better to consider, in the place of any one of the principal axes in particular, 
the instantaneous axis, which is a line inclined to these at angles a, B, y, the cosines of 


which are 2, 2, 7 (if as before w?=p?+q?+7%). Considering the instantaneous axis 


as a point of the sphere, let j denote the distance OP from the fixed point P, and 
$ the inclination OPQ of this distance to the fixed are PQ. We have 


cos 7 = cos a cos l + cos 8 cos m + cos y cos m, 
sin j cos $ = cos a sin l cos A+ cos £ sin m cos p + COS y Sin n COS v, 


sin j sin $ = cos a sin / sin X + cos f sin m sin w+ cos y sin n sin v, 


cos a — cos } cos) a cos y cos m — cos B cos n 
OMNE MT aie P NP MET E t 
cos 8 — cos m cos j COS a COS n — COS y cos | 
cos rr = T WES Spies Rua cos añ ER Te E A IÓN 
stills, sinmsinj  ” ram sin m sin j í 
COS y — cos n cos 7 cos 8 cos | — cos a cos m 
dos (y — y AAA Er DREEET HL Om. 


Sın n SIM ) sia n SIN ) 
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so that X, p, v are determined in terms of j and ¢. These expressions give 


d$ = Fan ig . {cos L (qdr — rdq) + cos m (rdp — pdr) + cos n (pdq — qdp)], 


which is reducible to Eulers equation 


p (MG — NB) cos l +q (NA — LG) cos m+ r (L38 — 174) 0058; 


d$ =dt E-23LMNu-—wv 


and thence, substituting for cosl, cosm, cos» their values, and observing that 
Ap (MG -— NL) + B? (N9IL- LG) + C? (DB — MA) 2 (H - 2LMN Fu), 
BCp (MG — NBF + CAP (N9L — LG? + ABr (LB — MAY = F (H —2LMN Fu), 
LA (MG -JN38) + MB (N3L— - LG) +NC (LY — MA) = LMNF, 
the equation becomes 


- 9H - 2LMPu) , PH 2LMN Pu) y (G — 9») 
G GV (K—2LMNGu) 


, IMNFpqr v (G - D) D) 
Ts /(G) (K -2LMNGu) 


do (E —2LMNu — v’) + dt = sin U 


cos U, 


where it is to be remarked that 
G? (E—2LMNuw — v) 
=(G — D) F? + G (K —2LMN Gu) —(G —D*) (K — 2LMN Qu) sin? U 
— 2DF V (G — D’) y (K -2LMNGu) sin U. 
Now 


dt dy 2LMNFu)/(G) q, 


dU = —9LMNGu 


= pqrdt, 


the differential d$ can be expressed as a fraction, the numerator whereof is 
— DdU (K - 2LMNGu) y (G)+ Fd U y [G (G — D?) (K - 2LMNGu)} sin U 


y LMNFGdu y (G (G — D) 


UTN) 055 


and the denominator 
(G — D) F? + G(K - 2LMNGv) - 29 F y (G— D?) (K -2LMNGu)sin U 


—(G—D*) (K — 2LMN Qu) sin? U. 
To simplify, write 
y(K-2LMNGu)=s, /(G@—D)=h, 
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the numerator is 

— Ss dU y (G) -- Fhs dU y (G) sin U — Fhs y (G) cos U, 
and the denominator 

IF? + Gg —2DFhs sin U — ks sin? U, 
which, observing that h*=(G-— $*, is equal to 

(Fh —Ds sin Uy + Gs? cos? U, 


and we have 
— OsdU + Fhs sin UdU — Fhds cos U 
db = “Ci, — Deain OF + Gstcost O (5 


the integral of which is 
Fh —Dssin U 
ran vx aa! A ca FE J 
nee scos U y (G) ' 
where § is the constant of integration, or substituting for h, s their values, the 
equation is 
_ Fy (G&—D)-— D sin U y (K-2LMNGu) 
Gate tU = IA IN Ts 
It may be added that 


eibi ye H [DF + y (G — D) (K — 2LMNGu)] sin U), 


and therefore 
.. SF X ((G— D) (K — 2LMNGu)) sin U 
jp >  GX(E-2LMNu) 


Euler remarks that the complexity of the solution owing to the circumstance that 
the fixed point P is left arbitrary; and that the formule may be simplified by taking 
this point so that G — $?— 0, and he gives the far more simple formule corresponding 
to this assumption; this is in fact taking the point P in the direction of the normal 
to the invariable plane, and the resulting formule are identical with the ordinary 
formule for the solution of the problem. The term invariable plane is not used by 
Euler, and seems to have first occurred in Lagrange's * Essai sur le probléme de trois 
corps,” Pria de l'Acad. de Berlin, t. 1x., 1772. 


To prove the before-mentioned equation for d$; starting from the equations 


à Y 
cos j = cos a cos l + cos f cos m + cos y cosm = z, 


sin j cos $ =.cos a sin l cos A + cos £ sin m cos u + cos y sin n cos v 


sin j sin $ = cos a sin l cos A + cos B sin m sin y + sin y sin n sin n, 
we have 


cos j dj cos $ — sin j sin $ d$ 


= — sin a da sin l cos A — &c. + cos a cos X cos l dl + &c. — cos a sin ¿sin A da + &o., 
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the second term is 


P. cos X cot l (q cos n — r cos m) 
+ 2 cos y cot m (r cos l —pcosn) 


+ cosy cot n (p cos m — q cosl), 


and the third term is 
E p 


^ sin A cosec / (qcosm 4-rcosm) 
w 


+2 sin p cosec m (r cos n. + p cos l) 


w 


3 - sin v cosec n (p cos | + q cos m). 


Hence the second and third terms together are 


cos l cos n cosmcosn  . c my . cos l 
— I (cos p E — cos p PIRE sin X ——- + sin = ) + &e, 
e sin J sin m : sin 4 sin m 


Ei end ue — X) + sin X sin z sin (v — X) 
@ 


7 3 F 3 &o., 
+ cos p sin n cos (p — v).+ sin p sin n sin (u — E is 


2 dns k cos A cos (v — A) + sin X sin (v TH + &c., 


w + cos p cos (u — v) + sin usin (p, — v) 


ED PY sin n m cos bs dun = 2d + &o., 


P + cos [u — (u — v)) 


it 


sin n (— cos v + cos v) + &c., = 0; 
w 


we have therefore 
cos j dj cos $ — sin j sin $ d$ 


= — sin a da sin l cos X — sin B df sin m cos p — sin y dy sin n cos v,. 


— d P sinl cos A +d! . sin mcos y +d2. sin nsin v 
@ @ @ 
ya. 3 " 
cb (sin / cos A dp + sin m cos u dq + sin n cos v dr) 
do ,. 1 : : 
b (sin / cos X p + sin m cos y q + sin n cos v T) 


— — cot j cos $ d ^ — sinj sin $ d$. 
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Hence therefore 


sin j sin $d$ = — cot j cos gd — 
-i (sin / cos X dp + sin m cos p dq + sin n cos v dr) 
do ,. 1 : : 
p (sin / cos X p + sin m cos u q + sin n cos v r) 
= — cot j cos $ È (cos l dp + cos m dq + cos n dr) 
; dw 
+ cot j cos $, (p cosl + q cos m + r cos n) 
-i (sin / cos A dp + sin m cos y dq + sin n cos v dr) 


do ,. , : 
ers (sin / cos A. p + sin m cos p . q + sin n cosv . T) 


| 


Tu dE ((— cot j cos $ cos } — sin ¿cos X) dp + &c.} 


+ a {( cot ¿cos $ cos [L-- sinl cosX) p. + &c.]. 


But we have 


conu Wee Pd 


sinIsinj  ' 
= 5, — cot 1 cot j 
sin / sinj 


COS y cos m — cos B cos n 
sin ¿sin j 


sin (A — $) = 


3 


and thence 


OUS, cos (e. d) ddl cos $ (cos a — cos } cos j) — sin $ (cos y cos m — cos £ cos n) 


sin /sinj 
whence also 


cot j cos $ cos ¿+ sin cos X 


1 : i 
= -—.. [cos ¢ cos [ + cos $ (cos a — cos } cos j) — sin $ (cos y cos m — cos 8 cos n)}, 


sinj 


1 : 
= "E {cos a cos $ — sin $ (cos y cos m — cos B cos n)}, 


ji 


PT (p cos $ — sin $ (r cos m — q cos n)}. 
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Hence the expression for sin ¿sin $ de is 


1 
> fiir [(p cos $ — sin $ (r cos m — q cos n)) dp + ...] 


ue 


ee [{p cos @ — sin $ (r cos m — q cos n)) p+...] 


nu . [e dw cos $ — sin $ ((r cos m — q cos n) dp + ...]] 
ey Aue 
+ ans w? cos $ — sin j sin $ de, 
or finally 


sin j sin $ d$ = "- Tr [(r cos m — q cos n) dp + &c.], 


that is 


Hp 1 5) (r cos m — q cos «rie 


o^ sin’) 4 (p cos n—r cos 1) dq 
\+(q cos l — p cos m) dr j 


which is the required expression for d$. 


Recapitulating, 4, B, C, p, q, r denote as usual, 


B-—C C—A A—B 
bui , M= wt e Erant du — pqrdt, 
p —X (A + 2Lu), 
q =v (B +2Mu), 


r= (E +2Nu); 


A+, B+ G =£, 
NA +BB+C6C =F, 
9LA? + BB+ CC? = 


L-ABE + MBCA + MOAL =H, 
IABE + MBCA + NAB — KC ; 


so that 
K= EG- Fe, 
U=6+| (H —2LMNFu) du y (G) 
r ! (K — 2LMNGwu) y (A + 2Lu) (B + 2Mu) (€ + 2Nu)} ’ 
E YR 


19 
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cos U, 


_DAp , BOp(BG& — CB) (G — 3») . LAqr y (G — 3») 
esi =a" * GJ(K-3LMNGa  ** "+ eae -3LMNGw) 


_DBg CAg(CA— AG) (G — 9) . MBrp y (G — 3») 
wm= Cet sp ie Se U. 


DOr ABr(48 — BIDV(8—3*) a y, Npa y (C-D) 
G * Gy(K-3LMNGw) V(G(K - 2LMXNGu)) 


o?—E-—2LMNu, 


. SF ((G—99(K 2LMNGu) sin U 
Res m GEN Fa 


cos U, 


cos" = 


v = p cos l + q cos m 4- r cos 
= S [DF V ((G—D')(K - 2LMNGu) sin U), 


_ F4(G— 3») 9 sin Uy (K -2LMNGu) 
AA AL O 


[The angles which determine the position of the body are thus expressed in 
terms of u, which is given as a function of t by the foregoing equation du = pqrdt, 
where p, q, r denote given functions of u.] 


www.rcin.org.pl 


